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Radiative Transfer in Cometary Dust Atmospheres: Critique of Recent Developments 


Oldwig von Roos 
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Abstract* 

The gas and dust production rates of comets near the sun are intimately 
linked to the fate of the solar radiation as It Is scattered and partially 
absorbed by the cometary dust ••i.^elope or halo surrounding the comet's nucleus. 
This Is so, because the radiative energy Impinging on the surace of the nucleus 
governs the rate of evaporation of the nuclear matter (mostly water ice) and 
the release of dust particles embedded In the Ice which are entrained by the 
escaping gas and subsequently form a dust atmosphere or halo. On the other 
hand, the Interaction of the dust cloud with the solar radiation will In turn 
be responsible for the amount of radiative energy received by the nucleus and 
thus determine the evaporation rate. Here It will be shown, that recent attempts 
to quantitatively formulate this problem are In error. The correct formulation 
will be given and It will be pointed out that negative extinction (enhancement 
of the primary radiation received by the nucleus due to multiple scattering) 
may not be as large as predicted. 
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1.) Introduction. 

The Interaction of the solar radiation with a comet's nucleus as It 
approaches and recedes from perihelion, the bulld~up of a gaseous and dusty 
halo dwarfing the size of the nucleus, the Interplay between trapped radla*> 
tlon within the dust atmosphere surrounding the nucleus and the amount of 
radiative energy delivered to the nucleus determining Its rate of evapora- 
tion and production of dust by aerodynamic forces constitutes a topic of 
considerable research^). But It Is only recently that a comprehensive 
theory of radiative transfer within a come’:ary dust atmosphere has been 
promulgated^). Absorption and scattering of light by the dust particles of 
the coma are Important events, since they determine ultimately the amount 
of energy deposited on the surface of the cometary nucleus and thus, by 
this very fact, determine the density of the dust cloud and the gas produc- 
tion rate. In the chain of events, the radiative transfer throughout the 
dusty atmosphere of the comet constitutes one important link. 

Unfortunately, the theory given In ref. 2 contains a number of errors 
and it is the aim of this paper to correct these and to provide a general 
framework for the evaluation of the light-dust Interaction taking multiple 
scattering to arbitrary order Into account. There are essentially two ways 
at our disposal to deal with radiative transfer, l.e., the analysis of a 
radiation field in a medium which absorbs and scatters radiation. One 
way, the customary way^), consists of writing down the energy balance 
for the Intensity of the radiation (intensity i energy flux per frequency 
Interval and solid angle) thus: 



a I 
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+ e 
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Where ly la the Intensity of the radiation field at any given point In 
apace, a the path length over Which Che change of the Intenalty la followed, 
Oy the abaorptlon coefficient and ey the emlaalon coefficient and flnilly 
V la the frequency contemplated. In as much aa Cy dependa Itaelf on the 
scattered Intensity among other things, equation (1) constitutes a linear 
Integral equation. The solution of this Integral equation given appro** 
prlate boundary conditions can only be found approximately^). 

There exists however another method for dealing with radiative trans*' 
fer. We call it the multiple scattering method. It consists of following 
each ray as it enters the medium, is scattered once then twice and so on, 
summing up all contributions of all configurations in an infinite series 
and constructing the radiation field at each point in space in this manner. 
This method is of an advantage when either the absorption cross-sections 
for the radiation are large or the scattering cross-sections are small. 
It is also advantageous if the assumption of local thermal equilibrium 
cannot be made. It has been used recently for evaluating the interaction 
of solar radltlon with the cometary dust clouds^). While the first method 
(eq. (1)) leads to an Integral equation which Is uifflcult to solve, the 
second method leads to multiple Integrals which cannot be evaluated in 
closed form. 

An analogy exists between the two methods of dealing with radiative 
transfer and a problem in optics. The diffraction of light by a plane 
parallel plate and the evaluation of transmission and reflection coeffi- 
cients may be undertaken by either solving Maxwell's equations with the 
appropriate boundary conditions at the two surfaces of the plate or one 
follows a ray as it enters the plate via the first surface, being partially 
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reflected and tranenitted, the traneoitted ray being again partially reflected 
and transmitted by the second surface and so on. Adding up all Infinite many 
contributions leads to the same result as that obtained by the "boundary condition" 
method^). Obviously, the first method Is analogous to the Integral equation 
method (eq. (1)) and the second method to the multiple scattering method. 

In the next section we shall describe the multiple scattering method In 
detail and apply It to the Interaction of light with the dust cloud of a comet. 

2.) Analysis. 

Consider fig. 1. The cometary nucleus, a sphere of radius r)(, defines a 
polar coordinate system as shown. The z-axls or polar axis Is directed toward 
the sun. Therefore, the Incident flux of the solar radiation Eq In erg per cm** 
and sec Is pointed Into the negative z-dlrectlon. A dust cloud of number density 
N(r) per cm^ surrounds the nucleus and we assume that N Is only a function of r, 
the radial distance from the nucleus. For the formulation of the multiple scat- 
tering method this assumption Is not necessary but It simplifies the formalism. 
Ue assume with ref. 2 that Mle scattering by the dust particles (considered 

% 

spherical and of uniform density) constitutes the predominant mechanism by which 
the solar radiation Interacts with the dust cloud. Consequently, we Introduce an 
absorption cross-section o^(cm^l per dust particle and a differential cross- 
section for scattering dOg[cm^l defined by: 

db> 

dOg ■ o§ n (cos 6) An (2) 

where dot signifies the solid angle element Into which the Incident radiation Is 
scattered, 6 the angle between Incident and scattered radiation, the function n 
being a measure of anisotropy (for isotropic scattering n * 1 Independent of 6). 
Since 



A 



Where the Integration takes place over all solid angles, Og of eq. (2) 
signifies the total scatterii^ cross-section and therefore the extinction 
cross-section is given by: 


«ex - ®A + ®S» 

Of course, the solar flux and the cross-sections for Hie scattering (2) and 
and (4) are all very much frequency dependent. But since no inelastic 
scattering events take place we have refrained from explicitly showing this 
dependence; it is implicitly understood. Also, we assume, again for simpli- 
city, that the dust particles consist of one species only with constant 
radius a. The generalization to a di«t:ributlon of different grain sizes is 
trivial. 

Suppose the Intensity of radiation at position 1 and going into direc- 
tion 12 is I 2 (1). Then, the amount of intensity arriving at point 2 having 
not been scattered is given by: 

2 

1(2) - l 2 (l) expj- Ogx y* Ndzj , (5) 

1 

where the integration extends along the straight line between point 1 and 
point 2 and the integral is to be taken positive. 


In order to calculate the total amount of energy per second B received 
by the nucleus from the sun taking multiple scattering into account, we pro- 
ceed as follows. Without any dust cloud E would simply be: 

.it/2 


£ ■ £, 
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CO80 dto • 2itEg rjj / COS0 sln0 d0 • irrjj Eg, 


( 6 ) 


5 


since the angle between the nomal to a surface elenent of the nucleus and the 
direction of the Incoming radiation Is 6* However, with the cloud of dust rein- 
stated we may split the radiation Impinging on the surface of the nucleus Into 
parts of ever larger complexity* For Instance that part of the radiation which 
has undergone no scattering or absorption is given by (using eq. (5))t 


e 

E 

o 


E. 



COS0 exp 


I 



slnB dd d^ 


(7) 


We have put the upper limit of the integral over N to Infinity for convenience 
because N decreases rapidly with distance. From geometrical considerations It can 
easily be gathered that: (g) 


E - 2n 


r 

^ o 


dO cosO sln6 exp 


“ o 


ex 
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N((r^‘+ 


+ 2 r^ z 


CO80)^^^)dz ' 


and If the number density N Is kno%m as a function of Its argument, the Integral (8) 
Is also known In principle* For a quadratic dependence t>f N on r: 

N - n^/r^, (9) 


eq* (8) goes over Into: 
E 

where 


/ «/2 


COS0 sln0 d0 e 


-r0/sln0 


( 10 ) 


® " ®ex 


( 11 ) 


Expression (8) Is equivalent to eq* (5) of ref* 2* 
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Next we eonclder elngle ecetterlng evente. A volume element at poaitlon 
'l* ^1* ^1 magnitude d rj^ contains N(r^)d''r]^ duet particles* 
The incident solar radiation flux at this point has been attenuated and ia given 


hy: 


B(l) - Eg expj - 


N((r^^ + + 2rj x co80j)^^^)ds | <12) 


The radiative energy per sec ecattered into an arbitrary direction by the 
email amount of dust at this point is given by: 


dEj^ " Og ^ (co 8 «) ~ E(l) N(r^) d^r^. 


(13) 


where 5 is the angle between the incident direction and the direction of scat* 
taring. The amount of radiation impinging on the nucleus due to the radiation 
having cuffered one scattering event is given by eq* (13) when integrated over 
all possible positions jri ^<^d all directions from point Pj which intersect the 
nucleus (the cone emerging from P 3 in fig. 1 ) provided that properly taken into 
account is the attenuation of the scattered radiation as it emerges from point 
Pj^and strikes any part of the nuclear surface visible from P]^. It is now con- 
venient to introduce another spherical coordinate system with origin at Pj (P 3 in 
fig. 1 ) such that any line from. P]^ intersecting the nucleus makes an angle a 
with the line connecting the center of the nucleus with Pj which evidently 
became the z-axis of this new coordinate system. It is clear from fig. 1 (where 

P. replaces P. ) that 0 and ^ can take on all values between 0 and 6 and 0 and 

J i- max 

2 i> respectively, where 6 ^^ is determined by: 


sin 9 


max 




(lA) 


(r|(/r 3 in fig. 1). Furtharmore, th« aeattarlng angla 6 obeytt 


eosA ■ eo«9 eosS^ 4- sinO •in6|^ eos9 • 


( 15 ) 


The attenuation of the scattered radiation emerging from point P]^ and 
heading Into the direction 6, 9 In the new coordinate system Is given byt 


tj^cosd - (rjj^ - sln^e)^^^ 


I “ 6Xp V • 0 




^ - 2rj^ z co80)^^^) dz J • 


(16) 


Collecting eqe. (12), (13) and (16) and Integrating yields for that part 
of the energy Impinging on the nuclear surface which underwent one scattering 
event : 


f^\ 


N(rj) 




/ ^nax 

sine d0 J d(|) 


n(cos6). 


(17) 


According to fig. 1 no scattering event can take place In the shadow region, 
at least not for the first scattering event. Therefore, In the Integration over 
r^ the shadow region must be excluded. The prime Is a reminder of this fact. 
Analytically, this statement Is expressed by the definition: 

•2v\ ,it/2 


/ dV... /sme^ ‘‘’"r** ''r-* 

•' ^ •'o -'r V®!"® 

‘k K 


(18) 


For isotropic scattering (n ■ 1) eq. (18) goes over Into: 

• T ■ °ex/ "^^’^1^ ■*■ *^ ■*■ ^'l * co®®i>^^^>«*^ 


max 

slnee"^^®^de. 


(19) 


It Is eq. (19) Which must be compared with Its equivalent expression (16) of 
ref. 2. 
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If we set P(6) • 0 in the last (angular) integral of eq. (20) we obtaint 

'*®max 

sine de - 1 - yi-rj^ 2 /rj ^2 ^ (20) 




o 

and eq. (19) goes over into eq. (6) of ref. 2. But putting P(6) ■ 0 is equiva- 
lent to setting 0 ^^ ■ 0 according to eq. (16) Which means that the scattered 
light is not attenuated as it traverses the dust cloud toward the surface of 
the nucleus! When the dust particle density N is given by eq. (9), the inte- 
gral of eq. (16) can be performed and yields: 


F<6) - e(r^/r^sln6)(ir/2-e-tan”l^r^2/ri2*inViy/2)^ ^21) 
with e given by eq. (11). 

We now turn to higher order scattering contributions. Consider the light 
being scattered at Pi into the direction of ?2 (fig* D* The aonunt of radia- 
tion scattered into this direction is given by (see eq. (13)): 


where 


dE^ - Og n(cos6i) Ed) N(rj) d^r^^ , 
cosdj^ • ti2”^(ri cos6| - t2 cos02)* 


( 22 ) 

(23) 


The latter equation follows from simple geometrical considerations. ri2 is the 
length of the vector X2 ” Xl* 


with 


ri2 - (rj^^ + r2^ - 2rjr2 cosY22)^^^» 
cosYi 2 • cosOj^ cos02 ainOi sin02 cos(0i-02)* 


(24) 


( 25 ) 
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The avount of radiative flex arriving at ?2 will be attenuated in going 
fro® ?! to ?2 and ia given by: 



n(coedj^) E(l)N(rj^) d”'rj^ (4vrj^2^)“^ 


exp {-0 


./■ 


dz 


!■ 


where dA « rj^2^*l'*** presence of the factor rj^2**^ accounts for the divergence 
of the hean emerging from the point P^* At point ?2 the second scattering event 
takes place which may either lead to a third scattering event elsewhere (at P3 
in fig. 1) or scatter the light directly to the surface of the nucleus. The ana- 
lysis from this point on is identical to the one which lead to eq. (17), so we 
will not expatiate on it any further. The result is: 


(27) 




2 N(rj^)N(r 2 )K(l)n (cosdj^) exp 


r 2 I r®«ax 

-o^lN(r)d*J *^12 / ®l"®d8/dpn(co8 42^®^^*^' 


E(l) in (27) is given by eq. (12). determined by sinB^ax ” 

62 is the scattering angle between the direction r2 ^^d the B, ^ 
direction of the new polar coordinate system within the cone as discussed pre- 
viously is given by: 


COS62 ■ (tj cosY}^ 2 " *^2^*^12~^ ^*^1^*^12^ (slnB cos^ fslnB2 cosBj^ - cosB2 
slnBjcos(^2”^l)l slnB sln^ sinBj^ sin(^2~^l))* ^ 

The primes at the Integration signs indicate that the appropriate shadow zones 
must be excluded. For the Integration over rj it is shown in fig. 1 and for the 
integration over r2 it is shown in fig. 2. The analytic representation, simple in 
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the case of fig. 1 and given by eq* (IS), nay be formulated In the following 
manner in case of fig. 2: Given the ai'^gUs 63 and ^2t determine yx2 

eq. (25). If Y < */2 “ "Here Sgiax given by eq. (14), with replaced 
by r 2 , the integration over r 2 extends from r^ to infinity. If Y > */2 - Smax* 
the range of integration over r 2 extends from 


rK/sln(Y + 0niax) < r2 < «. (29) 


^1 ~ ®max * prescription (29) goes over Into eq. (18) as it must. 

Th& Integral (27) constitutes the second term in an infinite series of mul- 
tiple scattering contributions toward the insolation of the nucleus. It is by 
now clear how to generalize eq. (27) to triple scattering contributions for 
instance. We will refrain from doing so, since the expressions become quickly 
impractical for computer work and ref. 2 stopped at this point too. It is not 
very easy to compare our final expression (27) with the work of ref. 2, since 
the results there have not been put into closed form as we have done in eq. (27). 
However, eq. (26) and eq. (13) of ref. 2 should be equivalent. Comparing the 
two expressions, it becomes evident that indeed they are essentially equivalent, 
were it not for the fact that the divergence factor r ^^2 missing In ref. 2. 
But the radiation emerging from a point diverges; the flux becomes smaller the 
further away from the point source and this must be taken into account. As 
the points Pj and P 2 coalesce and rj2 * 0 ®q* (26) is not valid anymore, we 
are in the presence of the near field of the scattered radiation. Obviously, 
no contributions toward scattering have developed yet. We must exclude a small 
region about the singularity at r^^^ “ 0* Integration over is under- 
stood as a principal value as far as the singularity is concerned. For a qua- 
dratic dependence on distance of the dust density eq. (9), the attenuation 
along the path 1-2 is readily obtained and is given by: 


exp { - o 


ex 


f 


N(r)dz; » exp 


[-c(ri2 


r^/ri r^)(y^2^siny 


12 ^^ * 


(30) 
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Again G l8 defined by eq. (11). In ease of laot .opie scattering we obtain from 


eq. (27) the somewhat simpler fonmilax 

s 


witv owuawwtawfc ^ ^ 

h-°-tj »»> 


max 
slnB dS e 


-F(e) 


(31) 


3. ) Discussion. 

The first feu terms of a series of contributions to the total Insolation of 
a cometary nucleus surrounded by a dust cloud have been derived. Each subsequent 
contribution takes one more scattering event Into account than the previous one 
and the total sum determines the total Insolation* Of course, the terms of this 
series become very rapidly very complex and even high vpeed computers could not 
possibly cope with the higher order terms. The only hope for keeping computer 
times low and still obtain reasonable results. Is a rapid convergence of the 
series and this is assured when the absorption Is large or the scattering cross- 
section small. Just where we stand In this respect with the dust cloud may be 
seen from the following argument. Assume the dust cloud to obey the law (9). 
Assume also that there are approximately 20 dust particles per m^ at a distance 


of 1000 km from Che center of the nucleus 


This determines to be 2x10^^ cm~'. 


Assuming furthermore the dues to consist of spheres of Ollvln^) with a refractive 
index of slightly less than two and a radius of 1 pm, we obtain for the extinction 
cross-section^^ : 


5 —8 2 

Ogjj ■ ® ^ex " 6*3x10 ” cm‘ 


(32) 


For a nuclear radius ry of 3 km, we obtain for the quantity e 


t ■ ®ex "o/>^K ■ 0*06 (33) 

and we see that the attenuation as for Instance expressed in eq. (30) Is not very 
large. But It is easy to see that the general term Is proportional to c" and 
because of eq. (33) rapid convergence Is assured* 
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So far we have only considered the Irradiation of the nucleus by the sun* 
But* the nucleus Itself radiates too* Partially, through Its albedo, the nuclear 
surface will reradlate energy which will be subjected to oniltlple scattering In 
the dust atmosphere Just as the sun's radiation previously considered* Partially 
some of the solar energy will be used In melting the water lee and some of the 
energy will penetrate the Interior of the nucleus. Also, the nucleus will emit 
radiation by virtue of Its temperature. Because of the linearity of the radiative 
transfer equations, the nuclear radiation field due to reflection and temperature- 
radiation may be treated Independently from the solar radiation In the same manner 
as discussed In section 2 . One must only Insure that the proper energy balance Is 
maintained. Rut, this Is not the place tf discuss these Items any further. Our 
aim was to derive and describe a correct formulation of the radiative transfer 
theory based on the concept of multiple scattering events. We will therefore dis- 
cuss here briefly the two points In which the present theory differs from published 
results^). It Is clear from the exposition In section 2 , that the two points In 
question are: (1.) the neglect of attenuation of the radiation as It Is scattered 
toward the nucleus (the cone In fig. 1) and (2.) the disregard of the beam diver- 
gence, the flux of radiation In going from point Pj^ to ?2 for Instance (see fig. 1 ) 

2 

being Inversely proportional to the area rj^2 point ?2 where the second 
scattering event Is taking place. 

While the neglect of attenuation of the radiation during Its last leg of travel 
after having undergone n scattering events (points P2 and P3 of fig. 1 for In- 
stance) always Increases the amount of radiation the nucleus receives, the omission 
of the beam divergence tends also to overestimate the Insolation of the nucleus by 
scattered radiation If the extinction Is not too large since the volume sampled by 
the Integration over r2 Is large. For large absorptlvltles however, the Influence 
may be In the opposite direction because now smaller distance rj^2 weighted more 
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heavily on account of the decreaaa with distance of the attenuation factor 



(see eq» (31))* But In this case higher order contributions (double or more 
acatterlnga) tend to be small anyway. 

In conclusion, it must be said, that the numerical calculations and the 
numerical results given In ref. 2 are rather doubtful. If the extinction of 
solar radiation In the cometary dust cloud (both light scattering and absorp- 
tion) is of any significance to the gas and dust production and ultimately 
responsible for the net loss of matter during each perihelion passage of the 
comet, then, only new extensive numerical calculations with the aid of the 
correct theory, as propounded on the previous pages, will give an adequate 
picture of the situation. 
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Fig. 2. The integration over r, excludes the shadow zone ^ jlig 

lowet part of the cone with the tip at P, . 


